Relativistic models for quasielastic electron and neutrino-nucleus 
scattering 



Carlotta Giusti'' a and Andrea Meucci 1 



(N 

o 

O 

Q 

(N 



Dipartimento di Fisica, Universita degli Studi di Pavia and 
INFN, Sezione di Pavia, via Bassi 6 1-27100 Pavia, Italy 



Abstract. Relativistic models developed within the framework of the impulse approximation for quasielastic 
(QE) electron scattering and successfully tested in comparison with electron-scattering data have been extended 
to neutrino-nucleus scattering. Different descriptions of final-state interactions (FSI) in the inclusive scattering 
are compared. In the relativistic Green's function (RGF) model FSI are described consistently with the exclusive 
scattering using a complex optical potential. In the relativistic mean field (RMF) model FSI are described by 
the same RMF potential which gives the bound states. The results of the models are compared for electron 
and neutrino scattering and, for neutrino scattering, with the recently measured charged-current QE (CCQE) 
MiniBooNE cross sections. 
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1 Introduction 

The analysis of the recent neutrino-nucleus scattering cross 
sections measured by the MiniBooNE Collaboration II 1 121 
requires reliable theoretical models where all nuclear ef- 
fects are well under control. The energy region explored 
requires a relativistic description of the process, where not 
only relativistic kinematics is considered but also nuclear 
dynamics and current operators are described within a rel- 
ativistic framework. Within the QE kinematics domain, the 
treatment of the FSI between the ejected nucleon and the 
residual nucleus is an essential ingredient for the compari- 
son with data. 

Relativistic models for QE electron and neutrino-nucleus 
scattering, with different descriptions of FSI, are presented 
and compared in this contribution. 

In spite of many similar aspects, electron and neutrino- 
nucleus scattering are two different processes and it is not 
guaranteed that a model able to describe electron scatter- 
ing data will be able to describe neutrino scattering data 
with the same accuracy. The large amount of data available 
for electron scattering and all the theoretical and experi- 
mental work done over several decades, which provided 
a wealth of detailed information on nuclear structure and 
dynamics 0141 . make electron scattering the best avail- 
able guide to determine the predictive power of a nuclear 
model. A good description of electron scattering data can 
therefore be considered a necessary prerequisite for any 
model aimed to describe neutrino-nucleus cross sections. 

In the QE region the nuclear response is dominated by 
the mechanism of one-nucleon knockout, where the probe 
interacts with a quasifree nucleon which is emitted from 
the nucleus with a direct one-step mechanism and the re- 
maining nucleons are spectators. In the exclusive (e, e'p) 
reaction the outgoing nucleon is detected in coincidence 
with the scattered electron, the residual nucleus is left in 
a specific discrete eigenstate, and the final state is com- 
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pletely specified. In the inclusive (e, e') scattering only the 
scattered electron is detected, the final nuclear state is not 
determined, and the cross section includes all the avail- 
able final nuclear states. For an incident neutrino, neutral- 
current (NC) and charged-current (CC) scattering can be 
considered. In NC scattering only the emitted nucleon can 
be detected and the process is inclusive in the lepton sector 
(the final neutrino is not detected) and semi-inclusive in the 
hadron sector (the emitted nucleon is detected but the final 
nuclear state is not completely determined). In CC scatter- 
ing the inclusive scattering where only the final lepton is 
detected can be treated with the same models used for the 
inclusive (<?, <?') reaction. 

For all these processes the cross section is obtained in 
the one-boson exchange approximation from the contrac- 
tion between the lepton tensor, which, under the assump- 
tion of the plane-wave approximation for the initial and 
the final lepton wave functions, depends only on the lep- 
ton kinematics, and the hadron tensor, whose components 
are given by bilinear products of the matrix elements of the 
nuclear current between the initial and final nuclear states. 
Different but consistent models to calculate the hadron ten- 
sor in the exclusive and in the inclusive scattering have 
been developed and are outlined in the following sections. 



2 The exclusive (e, e'p) reaction 

The theoretical framework for the description of the ex- 
clusive (e, e'p) reaction is the distorted-wave impulse ap- 
proximation (DWIA). The model is based on the following 
assumptions M3I4I : 

i) An exclusive process is considered, where the residual 
nucleus is left in a discrete eigenstate « of its Hamiltonian. 

ii) The final nuclear state is projected onto the channel sub- 
space spanned by the vectors corresponding to a nucleon, 
at a given position, and the residual nucleus in the state n. 

iii) The (one-body) nuclear-current operator does not con- 
nect different channel subspaces and also the initial state 
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is projected onto the selected channel subspace. This is the 
assumption of the direct-knockout mechanism and of the 
IA. 

The matrix elements of the nuclear current are then ob- 
tained in a one-body representation as 



I f(q) I fn) 



(i) 



where /" the one-body nuclear current, is the single- 
particle (s.p.) scattering state of the emitted nucleon, <p„ the 
overlap between the ground state of the target and the final 
state n, i.e., a s.p. bound state, and the spectroscopic factor 
A n is the norm of the overlap function, that gives the prob- 
ability of removing a nucleon from the target leaving the 
residual nucleus in the state n. The s.p. bound and scatter- 
ing states are eigenfunctions of a non Hermitian energy- 
dependent Feshbach-type optical potential and of its Her- 
mitian conjugate at different energies. In standard DWIA 
calculations phenomenological ingredients are usually em- 
ployed: the scattering states are eigenfunctions of a phe- 
nomenological optical potential determined through a fit to 
elastic nucleon-nucleus scattering data and the s.p. bound 
states are usually obtained from mean-field potentials. 

The model can be formulated in a similar way within 
nonrelativistic DWIA [5| and relativistic RDWIA frame- 
works |6]. In RDWIA, calculations are performed with a 
relativistic nuclear-current operator and four- vector rela- 
tivistic wave functions for the s.p. bound and scattering 
states. Both the DWIA and the RDWIA have been quite 
successful in describing (e, e'p) data in a wide range of nu- 
clei and in different kinematics 14 6 17181 . 



3 Inclusive lepton-nucleus scattering 

In the inclusive scattering only the outgoing lepton is de- 
tected, the final nuclear state is not determined, and all 
elastic and inelastic channels contribute. This requires a 
different treatment of FSI with respect to the exclusive pro- 
cess, where FSI are described by a complex optical poten- 
tial whose imaginary part gives an absorption that reduces 
the calculated cross sections. The imaginary part accounts 
for the fact that in the elastic scattering, if other channels 
are open besides the elastic one, part of the incident flux is 
lost in the elastically scattered beam and goes to the other 
(inelastic) channels which are open. In the exclusive scat- 
tering only the channel n is considered and it is correct 
to account for the flux lost in the considered channel. In 
the inclusive scattering all the channels contribute, the flux 
lost in a channel must be recovered in the other channels, 
and in the sum over all the channels the flux can be re- 
distributed but must be conserved. If the inclusive cross 
section is obtained from the sum of all the integrated ex- 
clusive one-nucleon knockout processes, due to the inter- 
action of the probe with all the individual nucleons of the 
nucleus, and FSI are described by a complex optical po- 
tential with an absorptive imaginary part, the flux is not 
conserved. The use of a complex optical potential seemed 
inconsistent with the requirement of flux conservation and 
in many calculations based on the relativistic IA real po- 
tentials have been adopted to describe FSI. 

In the simplest approach FSI are simply neglected and 
the plane-wave approximation is assumed for the scatter- 
ing wave functions (RPWIA). In a different approach an 



optical potential is used but only its real part is retained 
(rROP). The rROP conserves the flux but, independently 
of its numerical results, it is conceptually wrong because 
the optical potential has to be complex owing to the pres- 
ence of open inelastic channels. The energy dependence of 
the optical potential reflects the different contributions of 
the inelastic channels which are open at each energy and, 
under such conditions, dispersion relations dictate that the 
optical potential must have a non vanishing imaginary part. 
In a different approach the scattering states are given by 
the same relativistic mean field potential considered in de- 
scribing the initial nucleon states (RMF) 091 101 . The RMF 
model gives a consistent description of bound and scatter- 
ing states, it fulfills the dispersion relations [11 j and also 
the continuity equation. 

In the Green's function (GF) model H 1 21 1 31 141 1 51 1 61 1 71 1 81 
FSI are described in the inclusive scattering by the same 
complex optical potential as in the exclusive scattering, but 
the imaginary part is used in the two cases in a different 
way and in the inclusive process it redistributes the flux in 
all the channels and the total flux is conserved. 

In the GF model, under suitable approximations, which 
are basically related to the IA, the components of the hadron 
tensor are written in terms of the s.p. optical model Green's 
function. This result has been derived by arguments based 
on the multiple scattering theory [11] or by means of pro- 
jection operators techniques within nonrelativistic [12 | and 
relativistic H 1 31 1 41 151 frameworks. The explicit calculation 
of the s.p. Green's function can be avoided 012113114 1 by 
its spectral representation, which is based on a biorthogo- 
nal expansion in terms of a non Hermitian optical-model 
Hamiltonian 'H and of its Hermitian conjugate 1-1' . The 
components of the hadron tensor are then obtained in the 
form|T3l 
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where Ef denotes the energy of the final nuclear state, s n 
denotes the eigenvalue of the discrete eigenstate n of the 
residual nucleus, and 



rf(£,£) = (<p n I j^(q) yll- f V'(E)\x { :\E)) 
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The factor Vl - "V'(E), where <V'(E) is the energy deriva- 
tive of the optical potential, accounts for interference ef- 
fects between different channels and justifies the replace- 
ment in the calculations of the Feshbach optical potential 
*V, for which neither microscopic nor empirical calcula- 
tions are available, by the local phenomenological optical 
potential 112I13L Disregarding the square root correction, 
the second matrix element in Eq. (01 is the transition am- 
plitude of single-nucleon knockout of Eq. ([1), where the 
imaginary part of the optical potential accounts for the flux 
lost in the channel n towards the channels different from 
n. In the inclusive response this loss is compensated by 
a corresponding gain of flux due to the flux lost, towards 
the channel n, by the other final states asymptotically orig- 
inated by the channels different from n. This compensa- 
tion is performed by the first matrix element in the right 
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hand side of Eq. ©, that is similar to the matrix element 
of Eq. ([TJ but involves the eigenfunction x^ E \E) of the 
Hermitian conjugate optical potential, where the imaginary 
part has an opposite sign and has the effect of increasing 
the strength. As a consequence, in the GF model the imag- 
inary part of the optical potential redistributes the flux lost 
in a channel in the other channels, and in the sum over n 
the total flux is conserved. 

The hadron tensor in Eq. © is the sum of two terms. 
The calculation of the second term requires the integration 
over all the eigenfunctions of the continuum spectrum of 
the optical potential. If the imaginary part of the optical 
potential is neglected, the second term in Eq. (fJJ vanishes 
and, but for the square root factor, the first term gives the 
rROP approach. 

In the usual applications of the GF model the matrix el- 
ements in Eq. (f2j) are calculated using the same phenomeno- 
logical bound and scattering states already adopted in DWIA 
and RDWI A calculations for exclusive one-nucleon knock- 
out reactions. For the sum over n a pure shell model de- 
scription is usually assumed: <p n are one-hole states in the 
target nucleus with a unitary spectra strength and the sum 
is over all the occupied states in the shell model. With this 
simplifying assumption the contribution of all the nucleons 
is correctly included in the inclusive response. 

The GF model allows to recover the contribution of 
non-elastic channels starting from the complex optical po- 
tential that describes elastic nucleon-nucleus scattering data. 
It provides a consistent treatment of FSI in the exclusive 
and in the inclusive scattering and gives also a good de- 
scription of (e.e') data 0121131191 . 

Both nonrelativistic GF H12I161 and relativistic RGF JT3] 
models have been considered for the inclusive electron scat- 
tering, while only the RGF has been used for CC neutrino 
scattering [14|. The results of the RMF and RGF models 
have been compared in 1 19 1 for the inclusive electron scat- 
tering and in lEOl for the inclusive CC neutrino scattering. 

An example is displayed in fig.[T] where the RGF, RMF, 
rROP, and RPWIA cross sections of the 12 C(e, e') reaction 
calculated in a kinematics with a fixed value of the incident 
electron energy (e = 1 GeV) and two values of the momen- 
tum transfer (q = 500 and 1000 MeV/c) are compared. 
Two parameterizations of the ROP have been used for the 
RGF calculations: the energy-dependent and A-dependent 
EDAD1 (RGF-EDAD1) and EDAD2 (RGF-EDAD2) flTH . 

The RGF and RMF results are always different from 
the results of the simpler RPWIA and rROP approaches. 
The differences between RMF and RGF, as well as the 
differences between RGF-EDAD 1 and RGF-EDAD2, in- 
crease with the momentum transfer. At q = 500 MeV/c 
the three results are similar, both in magnitude and shape, 
at q = 1000 MeV/c, the shape of the RMF cross section 
shows an asymmetry, with a tail extending towards higher 
values of the energy transfer a>, which is essentially due to 
the strong energy-independent scalar and vector potentials 
present in the RMF model. The asymmetry towards higher 
lo is less significant but still visible for RGF-EDAD 1 and 
RGF-EDAD2, whose cross sections show a similar shape 
but a significant difference in magnitude. At q — 1000 MeV/c 
both the RGF-EDAD 1 and RGF-EDAD2 cross sections are 
higher than the RMF one in the maximum region, but a 
stronger enhancement is obtained with RGF-EDAD 1, which 
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Fig. 1. Differential cross section of the 12 C(e, e') reaction for an 
incident electron energy e = 1 GeV and q = 500 and 1000 
MeV/c as a function of the energy transfer a>. Results for RPWIA 
(dotted), rROP (dot-dashed), RGF-EDAD 1 (solid), RGF-EDAD2 
(long dot-dashed), and RMF (dashed) are compared. 



at the peak overshoots the RMF cross section up to 40% 
and it is even higher than the RPWIA result. 

The differences between the RGF-EDAD 1 and RGF- 
EDAD2 results are basically due to the differences in the 
imaginary part of the two optical potentials. The real terms 
are very similar for different parametrizations and the rROP 
cross sections are not sensitive to the parameterization con- 
sidered. 

In fig. |2] the RGF, RMF, rROP, and RPWIA results are 
compared for the cross sections of the 12 C (v^,yLT) reac- 
tion calculated with the same incident lepton energy and 
momentum transfer as for the (e, e') reaction of fig. [1] For 
the RGF model, the RGF-EDAD 1 results are compared 
with the results obtained with the energy-dependent but A- 
independent EDAI potential (RGF-EDAI). Also in fig. [2] 
the shape of the RMF cross section shows an asymmetry 
with a tail extending towards higher values of a> (corre- 
sponding to lower values of the kinetic energy of the outgo- 
ing muon Tf,). An asymmetric shape is shown also by the 
RGF cross sections, while no visible asymmetry is given 
by the RPWIA and rROP results. Also in this case the dif- 
ferences between the two RGF results are due to the differ- 
ent imaginary parts of the two ROP's. As already shown 
for the (e, e') reaction, the RGF yields in general a larger 
cross section than the RMF. 

The results in fig. |2]present some differences with re- 
spect to the corresponding (e, e') cross sections in fig.Q] In 
both cases the differences between the results of the differ- 
ent models are generally larger for increasing value of the 
momentum transfer. For neutrino scattering, however, such 
a behavior is less evident and clear. In particular, the RGF- 
EDAD 1 cross section at q — 1000 Mev/c does not show the 
strong enhancement in the region of the maximum shown 
in fig.[T] where the RGF-EDAD 1 result is even larger than 
the RPWIA one. In the case of neutrino scattering the RGF 
results in the region of the maximum are generally larger 
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Fig. 2. Differential cross section of the ' 2 C(v,j,^~) reaction for E v 
= 1000 MeV and q = 500 MeV/c and 1000 MeV/c as a function of 
the kinetic energy of the outgoing muon T^. Results for RPWIA 
(dotted), rROP (dot-dashed) RGF-EDAD 1 (solid), RGF- EDAI 
(long dot-dashed), and RMF (dashed) are compared. 



than the RMF ones, but smaller than the RPWIA cross sec- 
tions. The numerical differences between the RGF results 
for electron and neutrino scattering can mainly be ascribed 
to the combined effects of the weak current, in particular 
its axial term, and the imaginary part of the ROP [20|. 

We recall that the RMF model uses as input the real, 
strong, energy-independent, relativistic mean field poten- 
tial that reproduces the saturation properties of nuclear mat- 
ter and of the ground state of the nuclei involved. As such, 
it includes only purely nucleonic contribution and does 
not incorporate any information from scattering reactions . 
In contrast, the RGF uses as input the complex energy- 
dependent relativistic optical potential. Phenomenological 
optical potentials, obtained through a fit of elastic proton- 
nucleus scattering data, are adopted in actual RGF calcu- 
lations. Therefore, the RGF model incorporates informa- 
tion from scattering reactions and takes into account not 
only direct one-nucleon emission, but all the allowed final 
states, as the s.p. Green's function contains the full propa- 
gator. 

The imaginary part of the ROP includes the overall ef- 
fect of all the inelastic channels, which give different con- 
tributions at different energies. This energy dependence 
makes the RGF results strongly dependent on kinematics. 
The differences between the RGF and RMF results can be 
ascribed to the inelastic contributions which are incorpo- 
rated in the RGF but not in the RMF (and in other mod- 
els based on the IA), such as, for instance, re-scattering 
processes of the nucleon in its way out of the nucleus, 
non-nucleonic A excitations, which may arise during nu- 
cleon propagation, as well as to some multinucleon pro- 
cesses. These contributions are not included explicitly in 
the RGF model with a microscopic approach, but they can 
be recovered, to some extent, by the imaginary part of the 
phenomenological optical potential. With the use of such 
a phenomenological ingredient, however, we cannot disen- 



tangle the role of different reaction processes and explain 
in detail the origin of the differences, but we can expect 
that the differences increase with the relevance of such in- 
elastic contributions. 

The comparison between the RGF and RMF results can 
therefore be useful to evaluate the relevance of inelastic 
contributions. If in many situations the two models give 
close predictions (usually different from those of the sim- 
pler RPWIA and rROP), there are also situations where the 
differences are large HI 91201 . 

In the comparison with data, we may expect that the 
RGF can give a better description of those experimental 
cross sections which receive significant contributions from 
non-nucleonic excitations and multi-nucleon processes. This 
is expected to be the case |22 23] of MiniBooNE data, 
given the nature of the experiment. While in electron-scattering 
experiments the beam energy is known and the cross sec- 
tions are given as a function of the energy transfer, in neu- 
trino experiments q and a> are not known and calculations 
for the comparison with data are carried out over the en- 
ergy range which is relevant for the neutrino flux. The flux- 
average procedure can include contributions from different 
kinematic regions where the neutrino flux has significant 
strength and processes other than direct one-nucleon emis- 
sion can be important [24 25 1. Part of these contributions 
are recovered in the RGF model by the imaginary part of 
the optical potential 



4 Comparison with CCQE MiniBooNE data 

The CCQE l2 C(v M , fi~) cross sections recently measured by 
the MiniBooNE collaboration [ 1 1 have raised debate over 
the role of the theoretical ingredients entering the descrip- 
tion of the reaction. The experimental cross section is un- 
derestimated by the Relativistic Fermi Gas (RFG) model, 
and by other more sophisticated models based on the IA, 
unless the nucleon axial mass M& is significantly enlarged 
(1.35 GeV/c 2 in the RFG) with respect to the accepted 
world average of all measured values ( 1 .03 GeV/c 2 0261271 ). 
mostly obtained from deuteron data. The larger axial mass 
obtained from the MiniBooNE data on carbon can also be 
interpreted as an effective way to include medium effects 
which are not taken into account by the RFG and by other 
models. Before drawing conclusions, it is therefore impor- 
tant to evaluate carefully the role played by all the nuclear 
effects. 

The effect of FSI is investigated in fig. [3] where the 
CCQE double-differential 12 C (v^,/0 cross sections av- 
eraged over the neutrino flux are displayed as a function 
of T/j for various bins of cos 9, where 6 is the muon scat- 
tering angle. The RMF results yield reasonable agreement 
with data for small angles and low muon energies, the dis- 
crepancy becoming larger as 9 and increase. The shape 
followed by the RMF and RGF cross sections fits well the 
slope shown by the data. The two models yield close pre- 
dictions at larger values of 7^, for all the bins of cos 6 
shown in the figure. The RGF cross sections are gener- 
ally larger than the RMF ones. The differences increase 
approaching the peak region, where the additional strength 
shown by the RGF produces cross sections in reasonable 
agreement with the data. The differences between the RGF- 
EDAI and RGF-EDAD 1 results are enhanced in the peak 
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Fig. 3. Flux-averaged double differential cross section per target 
nucleon for the CCQE ll C{v ll ,[T) reaction calculated in the RMF 
(blue line) and in the RGF with EDAD1 (red) and EDM (green) 
potentials and displayed versus T M for various bins of cos 6. The 
data are from MiniBooNE [ Q. The uncertainties do not include 
the overall normalization error &/V=10.7%. 



region and are in general of the order of the experimental 
errors. 

In fig.[4]the total CCQE cross sections per neutron ob- 
tained in the RMF, RGF, rROP, and RPWIA models are 
displayed as a function of the neutrino energy and com- 
pared with the "unfolded" experimental data [ 1 1 . The rROP, 
RPWIA, and RMF results usually underpredict the Mini- 
Boone cross section. It is shown in B28I29I that the differ- 
ences between models like the RMF, rROP, RPWIA, and 
superscaling tend to be washed out in the integration and 
that all these models, which represent essentially the same 
nucleonic contribution to the inclusive cross sections, un- 
depredict the total MiniBooNE CCQE cross section, whereas 
the inclusion of two-particle-two-hole meson-exchange cur- 
rents enhances the results. 

Larger cross sections, in particular for larger values of 
E Y , are obtained in the RGF with both optical potentials. 
The differences between the RGF-EDAI and the RGF-EDAD1 
results, being RGF-EDAI in good agreement with the shape 
and magnitude of the experimental cross section and RGF- 
EDAD1 above RMF but clearly below the data, are due to 
the different imaginary parts of the two ROP's, particularly 
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Fig. 4. Total CCQE cross section per target nucleon as a function 
of the neutrino energy E v calculated with the RGF-EDAD1 (solid 
line), the RGF-EDAI (long-dashed line), the RMF (dashed line), 
the rROP (dotted line), and the RPWIA (dot-dashed line). The 
experimental data for neutrino scattering are from MiniBooNE 
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for the energies considered in kinematics with the lowest 6 
and the largest T u . We notice that EDAI is a single-nucleus 
parameterization, which does have an edge in terms of bet- 
ter reproduction of the elastic proton- 12 C phenomenology 
irSTIl compared to ED AD 1 , and also leads to CCQE results 
in better agreement with data. 

The MiniBooNE collaboration has also measured 
CCQE events. The data analysis is currently ongoing [30 1 
and some preliminary results can be found on the Mini- 
BooNE website BTI . When available, the antineutrino mea- 
surements will be an additional source of information about 
the weak charged-current lepton-nucleus interaction. 

In fig- El the total CCQE cross sections per target nu- 
cleon calculated with the RGF, RPWIA, and rROP for an- 
tineutrino scattering are displayed as a function of the an- 
tineutrino energy E v . Also for v scattering the RGF re- 
sults are usually larger than the RPWIA and rROP ones. 
The differences between the RGF-EDAI and RGF-EDAI 
results are significant also for antineutrino scattering, al- 
though somewhat smaller than for neutrino scattering. More- 
over, we note that the antineutrino cross section does not 
saturate in the energy range up to w 2 GeV which we have 
considered. The different behavior of the cross sections 
calculated for neutrino and antineutrino scattering is re- 
lated to the relative strength of the vector-axial response, 
which is constructive in v scattering and destructive in v 
scattering with respect to the longitudinal and transverse 
ones 11321 . 

The MiniBooNE Collaboration has recently reported 
(2l also a measurement of the flux-averaged differential 
cross section as a function of the four-momentum trans- 
ferred squared, Q 2 = -q^q^, for neutral-current elastic 
(NCE) neutrino scattering on CH2 in a Q 2 range up to 
~ 1.65 (GeV/c) 2 . A careful analysis of v-nucleus NCE 
reactions introduces additional complications, as the final 
neutrino cannot be measured in practice and a final hadron 
has to be detected: the corresponding cross sections are 
therefore semi-inclusive in the hadronic sector and inclu- 



EPJ Web of Conferences 




E v [GeV] 

Fig. 5. Total CCQE cross section per target nucleon as a function 
of the antineutrino energy E v calculated with the RGF-EDAD1 
(solid line), the RGF-EDAI (dashed line), the rROP (dotted lines), 
and the RPWIA (dot-dashed line). 



sive in the leptonic one. Different relativistic descriptions 
of FSI are presented and compared with the NCE Mini- 
BooNEdatain 1551. 



5 Summary and conclusions 

A deep understanding of the reaction mechanism of neutrino- 
nucleus cross sections is mandatory for the determination 
of neutrino oscillation parameters. Reliable theoretical mod- 
els are required where all nuclear effects are well under 
control. Within the QE kinematics domain, the treatment 
of FSI is an essential ingredient for the comparison with 
data. 

Models developed for QE electron scattering and suc- 
cessfully tested in comparison with electron-scattering data 
have been extended to neutrino-nucleus scattering. The re- 
sults of different relativistic models for the inclusive QE 
electron and neutrino-nucleus scattering, with different de- 
scriptions of FSI, have been compared in this contribution. 
In the relativistic plane-wave impulse approximation (RP- 
WIA) FSI are simply neglected. In other approaches FSI 
are included in RDWIA calculations where the final nu- 
cleon state is evaluated with real potentials, either retaining 
only the real part of the relativistic energy-dependent com- 
plex optical potential (rROP), or using the same relativis- 
tic mean field potential considered in describing the initial 
nucleon state (RMF). In the relativistic Green's function 
(RGF) model FSI are described in the inclusive scattering 
by the same complex optical potential as in the exclusive 
scattering, but the imaginary part is used in the two cases 
in a different way and in the inclusive process it is respon- 
sible for the redistribution of the flux in all the channels 
and the total flux is conserved. 

The differences between the results of the different mod- 
els depend on kinematics. The predictions of the RGF and 
RMF models are close in many situations, and generally 
different from the RPWIA and rROP results. There are, 
however, also situations where the differences are large. 
Larger cross section are generally obtained with the RGF 



model, which is able to give a better description of experi- 
mental data. In particular, the RGF model is able to repro- 
duce the CCQE MiniBoone cross sections without the need 
to increase the standard value of the nucleon axial mass. 
The enhancement of the RGF cross sections is due to the 
translation to the inclusive strength of the overall effect of 
inelastic channels, due, for instance, to re-scattering pro- 
cesses, non-nucleonic excitations, which may arise during 
nucleon propagation, or to some multinucleon processes. 
Such inelastic contribution, which are not incorporated in 
the RMF and in other models based on the impulse ap- 
proximation, are not included explicitly in the RGF model 
with a microscopic approach, but they can be recovered, 
to some extent, by the imaginary part of the optical poten- 
tial. The use of phenomenological optical potentials, how- 
ever, does not allow us to disentangle the role of different 
reaction processes and explain in detail the origin of the 
enhancement. 

Other models, where multinucleon components are ex- 
plicitly included, are also able to describe the MiniBooNE 
data without increasing the value of the axial mass [34 35 36 37 38 1. 
The important role of contributions other than direct one- 
nucleon emission is therefore confirmed by different and 
somewhat alternative models. A careful and consistent eval- 
uation of all nuclear effects is required before definite con- 
clusions can be drawn. A detailed comparison of the mod- 
els and of their results would be helpful for a deeper under- 
standing. Processes involving two-body currents, whose 
role has been discussed in [25 27 28 29 1, should also be 
taken into account explicitly and consistently in a model to 
clarify the role of multinucleon emission. Fully relativistic 
microscopic calculations of two-particle-two-hole (2p-2h) 
contributions are extremely difficult and may be bound to 
model-independent assumptions. 

The RGF results are also affected by uncertainties in 
the determination of the phenomenological optical poten- 
tial. At present, lacking a phenomenological optical po- 
tential which exactly fulfills the dispersion relations in the 
whole energy region of interest, the RGF prediction is not 
univocally determined from the elastic phenomenology. The 
differences of the results obtained with different parametriza- 
tions of the relativistic optical potential are produced by the 
different imaginary part, which is the crucial ingredient of 
RGF calculations. It is interesting to notice that the best 
predictions in comparison with data are given by the EDAI 
potential, that is also able to give the best description of the 
elastic proton- 12 C phenomenology. A better determination 
of a phenomenological relativistic optical potential, which 
closely fulfills the dispersion relations, would be anyhow 
desirable and deserves further investigation. 
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